A0-A04a  902 


UNCLASSIFIED 


national  bureau  of  standards  WASHINGTON  D C F/6  20/3 

determination  of  charge  or  polarization  distribution  across  POL— etc (U) 
SEP  77  A S DEREGGlt  C M GUTTHaN.  F 1 HOPSIK  N00014-7e-F>0015 


1 OF  1 

^>0404890? 

a 

1 “ 

i «-  — 

1 

1-  ^ 

't/T. 

END 

DATE 

FILMED 

2-78 

ODC 

&) 


OFFICE  j)F  NAVAL  RESEARCH 


.^Reggl,  CtM.^uttnan.  F.I./Mopslk/ 
] G.T.^avlSjijii  M.G.^roadhurst  i 

Prepared  for  Publication  In 

Physical  Review  Letters 


National  Bureau  of  Standards 
Bulk  Properties  Section 
^Washington,  D.C.  20234  y ^ 

(7/7  / 


Reproduction  In  whole  or  In  part  Is  permitted  for  any  purpose  of  the 
United  States  Govemnent. 


Approved  for  Public  Release;  Distribution  Unlimited. 


D D C 

PTOPimfTT 
JAN  ta  isn 

ElSEDin 

' B 


ft 


I 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  /wiiwi  DM*  Enl*r*iO 

REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTION.S 
BEFORE  COMPLKTINO  FOKM 

I.  REPORT  NUMBER  / j.  OOVT  ACCESSION  NO. 

Technical  Report  #6 

S.  RECIPIENT'S  CATALOG  NUMBER 

4.  title  (end  SuHlUo) 

"Determination  of  Charge  or  Polarization 
Distribution  Across  Polymer  Electrets  by  tlu 
Thermal  Pulse  Method  and  Fourier  Analysis  ' 

S.  TYPE  OF  REPORT  S PERIOD  COVE^  . 

Technical  Report  #6 

f.  PERFORMING  ORO.  REPORT  NUMBER 

7.  AUThORI*! 

A.S.  DeReggi,  C.M.  Guttman,  F.I.  Mopsik, 

G.T.  Davis*  and  M.G.  Broadhurst 

1 

S.  CONTRACT  OR  GRANT  NUMBERf.l 

N00014- 7^0015'  ^ 

S.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

National  Bureau  of  Standards 

Bulk' Properties  Section,  Polymers  Division 
Washington.  D.C.  20234 

to.  program  ELEMENT,  PROJECT,  TASr 
AREA  4 WORK  UNIT  NUMBERS 

NR- 356-522 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Office  of  Naval  Research 

Chemistry  Program 

Arllnaton.  Va.  22217 

It.  report  DATE 

September  1977  ( 

IS.  NUMBER  OF  PAGES 

13 

u.  MONITORING  AGENCY  NAME  f AOpRCSS^II  fl/l^fwil  ffpm  CcntroiUng  OMIf) 

IS.  SECURITY  CLASS,  (ol  IMo  toport) 

Unclassified 

l».  DISTRIBUTION  STATEMENT  (ol  Ull»  K*po€t) 

According  to  attached  distribution  list. 


W 

1S«.  DECLASSinCATION/OOWNCRADINa 

schedule 


TOSiHIBUnON  statementT 

Appeared  for  public  releosf] 
OlstxlbQllon  Unllmltad 


( 

•7.  DISTRIBUTION  STATEMENT  (ol  Iho  oboirmcl  MtMatf  in  Block  20,  II  OlllotonI  horn  Repott)  | 


IS.  SUPPLEMENTARY  notes 


To  be  published  in  Physical  Review  Letters 


If.  KEY  WOllOt  fCoatlmi#  «n  tflf*  II  Itfpfilliy  ff  ffocJr 

Charge  distribution;  Fourier  analysis;  polarization  distribution; 
polyvlnylidene  fluoride;  pyroelectricity;  themal  oulse;  transient 
response;  vinylldene  fluoride  copolymer 


ABSTRACT  iCotillnu*  mi  tovoroo  olOo  If  itotooooir  anB  lOonlllf  tf  klock  miaik»r> 

Fourier  analysis,  applied  to  the  recently  introduced  thermal  pulse  method, 
yields  new  and  unique  relations  between  the  time- dependent , pyroelectric 
response  o£..a^  thin  specimen  and  the  Fourier  coefficients  of  the  charge 
or  polariz'AeU^ ‘distribution  across  its  thickness.  The  new  analysis  is 
applied  illustratively  to  thermal  pulse  data  for  a vinylldene  fluoride 
copolymer  electret. 


Decenlnatloii  of  Charge  or  Polarlaatlon  Distribution 
Acroas  Polyner  Electrets  by  the 
Themal  Pulse  Method  and  Fourier  Analysis 

by 

A.S.  DeReggl,  C.M.  Guttiun,  F.I.  Mopslk,  G.T.  Davis  and  M.G.  Broadhurst 

National  Bureau  of  Standards 
Washington,  D.C.  20234 

Abstract 

Fourier  analysis,  applied  to  the  recently  Introduced  thermal  pulse 
method,  yields  new  and  unique  relations  between  the  time-dependent, 
pymulectrlc  response  of  a thin  specimen  and  the  Fourier  coefficients 
of  the  diarge  or  polarization  distribution  across  Its  thickness.  The 
new  analysis  Is  applied  Illustratively  to  thermal  pulse  data  for  a 
vlnylldenefluorlde  copolymer  electret. 


□ n 


Several  new  experimental  methods  have  been  proposed  In  recent 


(123)  ’ 

Letters  ’ * to  probe  the  spatial  distribution  of  charge  and/or  polariza-  j! 

tlon  across  the  thickness  of  thin  (thickness  of  order  10  urn)  dielectric  samples. 

These  methods  are  currently  of  interest  as  means  to  study  piezo-  and  pyro-  ;i 

electric  polymers.  The  electron-beam  method  of  Sessler  et  al.^^^  has  i 

a resolution  of  _1  pm,  but  is  destructive  In  the  sense  that  the  measuring  j 

process  irreversibly  changes  the  quantities  measured.  Acoustical  methods  of  j 

(2)  I 

generating  a propagating  pressure  step  across  a sample  in  order  to  measure 

its  inhomogeneous  piezoelectric  response  are  non-destructive;  but  the  step  rise 

time  requirements  become  severe  if  a resolution  of  1 pm  is  desired  (a  resolution  i 

of  1 pm  requires  a step  rise  time  of  S 0.5  ns  assuming  a sound  velocity  of 

2000  m/s.)  Thermal  methods  of  producing  transient  inhomogeneous  strain  across 

a sample  in  order  to  measure  its  inhomogeneous  pyroelectric  response  as,  for  example, 

in  the  Collins's  thermal  pulse  method^^*^'^^  are  both  non- destructive  and 

comparatively  easy  to  implement  experimentally  (thermal  equilibration  times 

••J 

are  ~1  ms).  However,  because  the  thermal  equilibration  process  is  diffusive 

rather  than  propagative  (wavelike),  information  contained  in  the  experimental 

data  about  charge  and/or  polarization  distributions  is  convoluted  with  a time-  p, 

dependent  temperature  distribution  and  somehow  mtist  be  deconvoluted.  The  major 

and  often  quoted^^’^’^*^’^^  objection  to  the  Collins  thermal  pulse  method  has 

been  that  the  deconvolution  procedure^^’^*^^  yields  non-unique  charge  and/or  I 

polarization  distributions. 


In  this  Letter,  we  show  that  the  Collins'  experiment  yields  the  first  few 
Fourier  coefficients  of  the  charge  or  polarization  distribution,  not  the  charge 
or  polarization  di retribution  Itself  as  may  be  supposed  from  Collins'  papers 
Further,  we  present  expressions  for  the  response  which  allow  the  Fourier  coefficients 
to  be  determined  from  the  experimental  data  by  numerical  or  graphical  analysis.  We 
apply  the  new  expressions  to  thermal  pulse  data  for  a copolvmer  of  vinylldene 
fluoride  with  tetrafluoroethylene.  This  example  Illustrates  both  the  power  and 
the  limitations  of  the  thermal  pulse  method. 

Consider,  first,  the  temperature  T(x,t)  at  time  t across  a samole  regarded  as  a 
slab  of  Infinite  extent  and  of  thickness  d (see  Inset  In  Fig.  1)  after  a radiant  hear 
pulse  of  duration  t^  has  been  applied  to  the  surface  x-0.  Assume  that  1)  the 
aluminum  electrodes  have  negligible  optical  transmission  and  negligible  thermal  mass; 
2)  heat  flows  In  the  x-dlrectlon  only;  and  3)  the  sample  does  not  lose  heat 
to  the  surroundings.  The  Initial  temperature  T(x,0)  taken  at  the  end  of 
the  thermal  pulse  Is  written 

T(x,0)  - Tj^  -f  AT  (x,0)  (1) 

where  T^  Is  the  uniform  temperature  of  the  sample  prior  to  the  pulse,  and 

AT(x,0)  Is  the  deviation  from  T^.  In  practice,  ATCx,0)  Is  a sharply  localized 

function  extending  from  x-0  with  a width  C « d.  The  temperature  at  t > 0 Is 

(6) 

T(x,t)  ” ”^1  AT(x,t)  where 


AT(x,t) 


ao  2 

a + El  a cos(nTrx/d>  exp(-n  t/x,)* 
o n*-*-  n 1 


(2) 
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where  • d /(ir  k)  where  k Is  the  thermal  dlffuslvlty,  and  where 
1 

a^  ■ I AT(x,0)  dx  ■ lim  AT(x,t) 

J o t-x* 

and 

2 

®n  “ d ^T(x,0)  co8(nxx/d)  dx,  forn-1,2,... 

' n 


In  the  limiting  case  where  AT(x,0)  6(x)  all  a ’s  are 


a « 2a  , for  n-1,2,... 
n o 

However,  for  a general  AICx,U),  the  temperature  at  the  surfaces 
are  described  by 

AT(0,t)  "a  + ?■,  s exp(-n^t/T-) 

’ o n-1  n !•' 

and 

AT(d,t)  - a^  + „?iC“l)”  exp(-n^t/Tj^)  . 


Expressions  (6)  and  (7)  have  been  used  previously  as  the  basis  for  thermal 
dlffuslvlty  measurements In  the  present  case,  the  dimensionless  quantities 
AT(0,t)/a^  and  AT(d,t)/a^  can  be  obtained  by  measuring  the  transient  resistance 

(A) 

of  one  or  both  electrodes  as  done  by  Collins  \ Then  the  a /a  's  and  t,  can  be 

n o 1 

determined  without-  knowing  the  detailed  shape  of  the  light  pulse. 


r 


- 4 - 


Consider,  now,  the  inhomogeneous  pyroelectric  response.  Assume  that 
real  charge  of  volume  density  p(x)  and  elementary  dipoles  responsible  for 
polarization  ?(x)  - ik  P(x)  move  reversibly  when  the  sample  thermally  expands 
or  contracts . Following  Collins , ^ ^ the  open—clrcult  voltage  due  to  charge , 
4V^,  and  that  due  to  polarization,  AV  , are 


o -o  fd 


o ■'  o 


[ p(x)  f*AT(x',t)  dx'  dx. 


a + a - a fd 

AV  (t)  - -2 5 L I 

P ee  I 

o 


( P(x)  AT(x,t)  dx, 
* 


where  o^,  a^,  and  are  temperature  coefficients  (1/x) (dx/dT) , (1/e)  (de/dT), 

and  (1/P)  (dP/dT)  respectively,  e Is  the  relative  permittivity  and  c Is  the 

o 

permittivity  of  vacuum.  Substitution  of  (2)  into  the  above  equations  gives 


where 


I-o  ^ ^ 7 nil  e*p(-n2t/tpi. 


A - X p(x)  dx 
' o 


A - p(x)  pln(n«x/d)  dx,  for  n*l,2,... 
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AV^(t) 

P 


a + a - a ^ 2 

—2 la  B + E,  a • B exp(-n  t/x.)], 

ee  o o n»l  n n l 


where 


B - P(x)  dx 


/■d 

B - P(x)  cos(mrx/d)  dx  for 

^ ' r\ 


It  should  now  be  noted  that  a , A and  B are  coefficients  of  Fourier  series 

n n n 

expansions  for  ^iT(x,0),  p(x)  and  P(x)  respectively  if  these  are  expanded  as 
cosine,  sine  and  cosine  series  respectively. 

The  foregoing  analysis  shows  that  the  Inhomogeneous  pyroelectric  response 

2 

consists  of  exponentially  decaying  components  with  relaxation  times  x /n  . 

^a  A 

Measurement  of  this  response  yields  Fourier  coefficient  combinations  (— ) (t^) 

on  a A 

d ^n  *n  o o 

(— ) In  the  charge  case,  or  (— ) (— ) in  the  polarization  case  (or  a linear 

o o 

combination  of  these  in  the  mixed  case)«  For  sufficiently  short  pulse  duration 
2 

(t  <<  x./n  up  to  the  highest  n of  Interest),  the  approximation  a /a  » 2 

no 

can  be  made.  Otherwise,  a^/a^  can  be  obtained  from  a separate  measurement  of 

the  transient  resistance  of  the  electrodes.  Thus,  Fourier  coefficients 

of  the  charge  or  polarization  distribution* A /A  or  B /B  , can  be  determined  up  to 

n o n o 

terms  with  n of  order  N ■ /xj/t^.  The  spatial  resolution  of  the  thermal  pulse 
experiment  Is  roughly  d/N. 
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Limiting  cases  of  charge  and  polarization  distributions  are  of  Interest. 

Note  that  if  p(x)  6(x),  then  AV^Ct)  -*-0  (l.e.,  charge  on  surface  x - 0 does 

not  contribute  to  AV  (t)),  while  if  p (x)  -►  6(x-d),  AV  (t)  becomes  Indlstlng- 

P p 

uishable  from  AV^Ct)  for  uniform  polarization.  This  is  a manifestation  of 

the  electrostatic  equivalence  of  p(x)  and  -dP/dx.  Note  also  that  if  P(x)  ->  6(x),  then 

AVp(t)  becomes  proportional  to  AT(Ott)  in  (7),  while  if  P(x)  -►  6(x-d),  then 

AVp(t)  becomes  proportional  to  AT(d,t)  in  (8).  This  is  a consequence  of  the 

symmetrical  roles  played  by  a and  B . 

n n 

Figure  1 shows  signals  proportional  to  open-circuit  voltages  for  two 

copolymer  samples  50  pm  thick  Cvlnylidenef luorlde  with  27%  bv  weight  of 

tetraf luoroethylene) . For  each  sample  Curve  A is  for  light  pulsed  at  x-0 

and  Curve  B is  for  light  pulsed  at  x~d.  The  difference  in  the  samples  is 

in  the  poling  treatment.  Both  samples  were  stored  for  long  periods  under 

short-circuit  conditions  after  poling.  Sample  (a)  was  poled  at  room  temper- 
ed) 

ature  by  corona  charging  ' with  a bias  voltage  of  -2000  V.  Curve  A,  shown 

in  Figure  1(a)  is  nearly  proportional  to  the  AT(0,t)  signal  (not  shown)  at 

times  t >>  t^where  t^  « 20  ps  while  Curve  B is  nearly  proportional  to  AT(d,t). 

These  observations  indicate  (assuming  that  all  B *0  for  n>N)  a sample  with  strong 

n 

localization  of  polarization  towards  x*0  as  shown  qualitatively  by  (a)  in  the 
inset  (or  a sample  with  a charge  layer  Imbedded  Just  inside  the  sample  close  to 
x-0) . Sample  (b)  was  poled  during  crystallization  by  maintaining  300  V across 
the  sample  as  it  was  cooled  from  130°C  and  crystallized  for  6 hours  at  118. 5°C 
before  cooling  further  to  100“C  (0.3  K/min)  followed  bv  quenching  to  room 
temperature.  The  signals  for  this  case,  shown  in  Figure  1(b),  indicate  a 
distributed  asymmetric  pyroelectric  activity. 


L 


I 

i 

J 


1 
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Note  that  series  (10)  and  (13)  acquire  alternating  signs  If  the 

direction  of  Irradiation  Is  changed  (In  which  case  a changes  to  (-l)“a_) 

n n 

or  If  the  sample  Is  reversed  (In  which  case  A and  B change  to  -(-1)°A 

n n n 

and  -(-1)"b^).  Thus,  If  separate  response  signals  are  obtained  by  thermally 
pulsing  opposite  surfaces  In  turn,  the  even  and  odd  Fourier  coefficients  can 
be  obtained  separately  from  the  sum  and  the  difference  of  the  two  signals. 

Visual  comparison  of  the  two  signals  also  reveals  whether  the  distribution 
of  pyroelectric  responslvlty  Is  symmetric  or  asymmetric. 

Figure  2 shows  semi-log  plots  of  the  curves  \ and  B from  Figure  1(b) 
together  with  (A+B)  containing  only  even  terms  of  the  response  and  (A-B) 
containing  only  odd  terms  of  the  response.  Assume  that 

the  responses  are  associated  with  polarization.  Graphical  analysis  of  the 

linear  part  of  the  odd  response  gives  ■ 4. 8 ms  (In  good  agreement  with  the 

value  expected  from  the  thermal  properties  and  the  thickness  of  the  film)  and 

2aj^Bj^  - 14.  Removal  of  this  linear  part  gives  a remainder  (Curve  D)  with  a 

slope  9/t^  Indicating  that  ^a^B^  =1.5  and  that  ^a^B^  • 0.  Similar  analysis 

of  the  even  response  gives  2a^B^  ■ 10.2,  ^a^B^  : 5 (from  Curve  C)  and  2a^B^  a 0. 

As  the  pulse  duration  of  t^  « 20  ys  should  allow  about  10  terms  to  be  observed 

(N  - /Xj^/t^  > 10) , and  terms  for  n 4 do  not  appear  In  the  pyroelectric  response 

(though  they  appear  In  the  x*0  resistance  signal).  It  Is  concluded  that  the 

for  n 4 are  negligible  compared  to  the  four  leading  coefficients.  The  values 

of  «2®2  then  corrected  to  the  values  2.675  and  0.775  respectively 

by  forcing  coefficients  a B , - a.B.,  a-B-,  - a-B_  when  Inserted  Into  (13) 

OO  .11  c j j 

to  yield  a curve  B'  In  Figure  2 within  experimental  uncertainty  of  measured 

Curve  B.  The  polarization  distribution  obtained  from  the  determined  coefficients 

assuming  a,  - a.  a.  - 2a  Is  the  distribution  (b)  In  the  Inset  of  Figure  1. 

X Z J o 
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The  samples  discussed  in  this  Letter  were  chosen  to  Illustrate  the 

theory.  Other  copolymer  samples  and  homopolymer  (polyvlnylldenefluorlde) 

samples  under  different  poling  conditions  can  eidilblt  nearly  uniform  polar- 
(9) 

Izatlon  . In  such  cases,  the  response  Is  nearly  step- like  and  B /B  ■ 0 

no 

for  all  n. 

The  present  work  Indicates  the  source  of  ambiguity  in  the  Collins'  decon- 
volution procedure.  The  thermal  pulse  data  (under  conditions  similar  to  those 
in  Collins'  experiments)  yields  no  more  than  10  or  15  coefficients  (based  on 
N - /xj^/t^) . The  Collins'  electrical  analog  sought  to  obtain  discrete  distri- 
butions characterized  by  20  adjustable  parameters.  Any  discrete  distribution  (of  which) 
there  are  many)  consistent  with  the  determinable  Fourier  coefficients  would 
reproduce  the  measured  transients  within  the  noise  In  the  data. 

Partial  support  of  this  work  by  the  Office  of  Naval  Research  is  gratefully 
acknowledged. 


FIGURE  CAPTIONS 


FIGURE  1 


Inhonogeneoiis  pyroelectric  response  of  two  copolymer 
samples:  (a)  sample  with  pyroelectric  activity  con- 
centrated towards  surface  x>0;  and  (b)  sample  with 
distributed  pyroelectric  activity.  Signals  marked 
A are  for  thermal  pulse  applied  to  surface  x-0  and 
those  marked  B for  thermal  pulse  applied  to  surface 
x-d.  Inset  shows  thermal  pulse  arrangment.  Polymer 
samples  have  vacuunh-deposlted,  reflecting  (partially 
absorbing,  non-transmlttlng)  electrodes.  Polarizations 
(a)  and  (b)  shown  In  the  inset  are  obtained  from  first 
few  Fourier  moments  of  corresoondlng  signals. 


FIGURE  2 


Graphical  analysis  of  response  of  Fig.  1(b).  A and  B, 

experimental  data;  A + B and  A - B are  even  and  odd 

response  (see  text).  The  relaxation  time  - 4.18  ms, 

coefficients  a^B^  - 5.1,  a^B^  - 7,  “ 2.675,  a^B^  - 

0.775,  SjjBjj  ■ 0 for  n ^ 4 are  determined  as  explained 

In  text.  B'  Is  a plot  of  Eq.(13)  with  a^  replaced  by 

(“l)**a  using  these  a_^B  's.  C and  D represent  the 
n n n 

n > 2 and  n - 3 terms  respectively. 


PI.'IURE  2 
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